Exclusive p° electroproduction on transversely polarized protons 
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The exclusive electroproduction of mesons was studied with the Hermes spectrometer at the 
Desy laboratory by scattering 27.6 GeV positron and electron beams off a transversely polarized 
hydrogen target. Spin density matrix elements for this process were determined from the measured 
production- and decay-angle distributions of the produced mesons. These matrix elements em- 
body information on helicity transfer and the validity of s-channel helicity conservation in the case 
of a transversely polarized target. From the spin density matrix elements, the leading-twist term 
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in the single-spin asymmetry was calculated separately for longitudinally and transversely polar- 
ized mesons. Neglecting s-cliannel lielicity changing matrix elements, results for the former can 
be compared to calculations based on generalized parton distributions, which are sensitive to the 
contribution of the total angular momentum of the quarks to the proton spin. 

PACS numbers: 13.60.Lc, 13.88.+e, 14.20.Dh, 14.40.Aq, 12.38.Qk 



Exclusive electroproduction of mesons can provide 
new information about the structure of the nucleon be- 
cause of its relation to generalized parton distributions 
(GPDs) [l|,[2,[l. In Ref. 0] it has been proven that the 
amplitude for hard exclusive electroproduction of mesons 
by longitudinal virtual photons can be factorized into a 
hard-scattering part and a soft part that depends on the 
structure of the nucleon and the produced meson. In the 
case of exclusive vector meson production, also the pro- 
duced meson is longitudinally polarized (in addition to 
the virtual photon being longitudinal). The amplitude 
for the soft part can be expressed in terms of GPDs. 

GPDs provide a three-dimensional representation of 
the structure of the nucleon at the partonic level, cor- 
relating the longitudinal momentum fraction of a parton 
with its transverse spatial coordinates. They are related 
to the standard parton distribution functions and nucleon 
form factors 0, [1, leading twist, meson produc- 

tion is described by four types of GPDs: Hi-9^ 
and E'''^ , where q stands for a quark flavour and g for a 
gluon. The GPDs are functions of t, x, and where t 
is the squared four-momentum transfer to the nucleon, x 
the average, and ^ half the difference of the longitudinal 
momentum fractions of the quark or gluon in the initial 
and final state. The quantum numbers of the produced 
meson determine the sensitivity to the various GPDs. In 
particular, at leading twist, production of vector mesons 
is sensitive only to the GPDs iJ « , , » , and Es . 




FIG. 1: The lepton scattering and hadron production 
planes together with the azimuthal angles and (jjs- 



The transverse target-spin asymmetry in exclusive 
electroproduction of longitudinally polarized vector 
mesons by longitudinal virtual photons is an important 
observable, because it depends almost linearly on the 
GPD E 0. This is in contrast to the unpolarized cross 



section, where the contribution of E is generally small 
compared to the contribution of H. At leading twist, the 
asymmetry is proportional to sm{(f)—(f>s), where (j) and 4>s 
are the azimuthal angles about the virtual-photon direc- 
tion of the hadron production plane and the transverse 
part St of the target spin, respectively, with respect to 
the lepton scattering plane (see Fig. [T]). 

The cross section and asymmetry for exclusive elec- 
troproduction e+p —> e' + p'^ +p' can conveniently be de- 
scribed using spin density matrix elements [^,[l3,[lH- By 
using the angular distribution of the produced vector me- 
son and of its decay products, as described by the polar 
and azimuthal angles -d and (p (see Fig. [2|), one can sep- 
arate the contributions of mesons with longitudinal and 
transverse polarization to the measured asymmetries. If 
s-channel helicity conservation (SCHC) holds, the lielic- 
ity of the virtual photon is transferred to the produced 
vector meson. In that case studying the asymmetry for 
the production of longitudinally polarized vector mesons 
is tantamount to selecting longitudinal virtual photons. 
Measurements have shown that SCHC holds reasonably 
well for exclusive electroproduction of p'^ mesons on an 
unpolarized target at Hermes kinematics • Thus in- 
formation on the GPD E can be obtained from measure- 
ments of the transverse target-spin asymmetry in exclu- 
sive p° electroproduction. 

Ultimately, these studies will help to understand the ori- 
gin of the nucleon spin, because it has been shown Q 
that the x-moment in the limit t — > of the sum of the 
GPDs H'^ and E'^ is related to the contribution J* of the 
total angular momentum of the quark with flavour q to 
the nucleon spin. 

In this paper, measurements of exclusive p*^ electropro- 
duction on transversely polarized protons are presented. 
For the flrst time, values of the spin density matrix ele- 
ments (SDMEs) and the transverse target-spin asymme- 
try for this process were determined. 

The data were collected with the Hermes spectrom- 
eter during the period 2002 - 2005. The 27.6 GeV 
Hera electron or positron beam at Desy scattered off 
a transversely polarized hydrogen target [14| of which 
the spin direction was reversed every 1 — 3 minutes. 
The average magnitude of the target polarization was 
\Pt\ = 0.724±0.059. The lepton beam was longitudinally 
polarized, the helicity being reversed periodically. The 
net polarization for the selected data was 0.095 ± 0.005, 
mainly because more data were taken with positive he- 
licity. 

Leptons were distinguished from hadrons with an av- 
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FIG. 2: The polar and azimuthal angles of the decay 7r+ of the 
p" in the rest frame. The positive z-axis is taken opposite 
to the direction of the residual proton, while the angle ip is 
defined with respect to the hadron production plane. 



erage efEciency of 98% and a hadron contamination of 
less than 1% by using the information from an elec- 
tromagnetic calorimeter, a transition-radiation detector, 
a preshower scintillation counter, and a Ring Imaging 
Cerenkov detector. Events were selected in which only 
one lepton and two oppositely charged hadrons were de- 
tected. 

In the event selection, the following kinematic con- 
straints were imposed: > 1 GeV^, > 4 GeV^, 
and —t' < 0.4 GeV^. Here —Q^ is the squared four- 
momentum of the exchanged virtual photon, W the in- 
variant mass of the virtual-photon proton system, and t' 
the reduced Mandelstam variable t' = t — to, where —to 
is the minimum value of —t for a given value of and 
the Bjorken variable xb- The average value of for 
the exclusive sample was 25 GeV^. The condition on 
t' was applied to reduce non-exclusive background. 

An exclusive event sample was selected by constraining 
the value of the variable 
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where Mx is the missing mass and M the proton mass. 
The measured Ai5 distribution, which includes con- 
straints on the invariant mass of the produced hadron 
pair as discussed below, is shown in Fig. [3l The peak 
around zero originates from the exclusive reaction. Ex- 
clusive events were selected by the requirement /S.E < 0.6 
GeV. This resulted in a total number of 7488 events. The 
background from non-exclusive processes in the exclusive 
region was estimated by using a Pythia6 Monte Carlo 
simulation [isl . in conjunction with a special set of 
Jetset fragmentation parameters, tuned to provide an 
accurate description of deep-inelastic hadron production 
in the Hermes kinematic domain 17, 3]. The simula- 
tion gave a very good description of the AE distribution 
in the non-exclusive region. The background fractions 
in the exclusive region varied between 7% and 23%, de- 
pending on the value of Q^, xb, or t' , with an average 
over all selected data of 11%. 




FIG. 3: The Ai5 distributions of the measured yield (num- 
ber of counts within the acceptance divided by the inte- 
grated luminosity) (dots) and a Monte Carlo simulation with 
Pythia6 of the non-exclusive background normalized to 
the same integrated luminosity (histogram). The kinematic 
cuts and the requirements 0.6 GeV < < 1.0 GeV and 

Mkk > 1.04 GeV were applied. The selected exclusive re- 
gion is indicated by the dashed area. 



The invariant mass of the two-hadron system M-^t^ 
was determined assuming that both hadrons are pions. 
Resonant tt+tt^ pairs, i.e., pairs produced in the decay 
p^ TT+TT^, were selected by the condition 0.6 GeV 
< A/tttt < 1.0 GeV. Contributions in the M^^ spectrum 
from the decay of a meson into two kaons were ex- 
cluded by requiring Mkk > 1-04 GeV, where Mkk is 
the invariant mass of the two-hadron system calculated 
assuming that both hadrons are kaons. After subtracting 
the simulated contribution from the non-exclusive tail in 
the region AE < 0.6 GeV and correcting for the non- 
constant acceptance with M^^^, the Mtttt spectrum for 
exclusive events was fitted with a p'^-peak plus a lin- 
ear background. For the shape of the p°-peak Soding 
and Ross-Stodolsky parametrizations were used. In both 
cases the resulting background was found to be negligible 
(0.7 ±0.5)%. 

In the analysis the recently developed formalism for 
electroproductionof a vector meson from a polarized nu- 
The cross section for exclusive p° 
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cleon was used 
leptoproduction is written as 



da 



1 



(27r)2 dxBdQ^dt 



dip d(j) dip d(cos dxB dQ^ dt 

^'^ W{xB,Q^,t,(t),<l)s,'P,^), 



(2) 



with Tp being a similar angle as 0s i but now defined 
around the direction of the lepton beam, and 



da 



dxB dQ^ dt 



.r.r-^^si^l (3) 



dar 
di 
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where r„ is the virtual photon flux factor in the Hand 
convention [l9j . e is the virtual-photon polarization pa- 
rameter, and dar/dt and duL/dt are the usual Xb,Q^, 
and t dependent 7*p cross sections for transverse and 
longitudinal virtual photons, respectively. 

The function W{xb,Q^, t, 4>, (j)s, f, "d) describes the an- 
gular distribution of both the produced p° and its decay 
pions. It consists of several terms corresponding to dif- 
ferent polarizations of the incoming lepton beam and the 
target nucleon: 



W = Wuu + P^Wlu 
PiSlWll 



SlWul + 

StWut + PiStWlt. (4) 



where the left (right) subscript specifies the beam (tar- 
get) polarization: unpolarized (C/), longitudinally (L), or 
transversely (T) polarized, and Pf,, Sl, and St represent 
the longitudinal polarization of the beam, and the lon- 
gitudinal and transverse polarization of the target (with 
respect to the virtual photon direction), respectively. 

For the case of zero beam polarization and only trans- 
verse target polarization^ the angular-distribution func- 
tion reads 

W(0, (/.s, v', ^) = Wuui<P, V. ^) + StWut[<P, c^s.^, ^) ■ 

(5) 

Here and in the following the dependence of the vari- 
ous angular distribution functions W on xb, , and t is 
omitted for the sake of legibility. 

The functions Wjjy (with Y = U,T) can be further 
decomposed into terms corresponding to specific po- 
larizations, indicated by the superscripts, according to 



V2 cos d sin d W^l (0s , 0, (y9) + sin^ d W^^ {4>s,(t>,v) 

(6) 

Note that in the case of Wjju there is no dependence 
on (jjs- The production of a longitudinally polarized 
is described by W^y^ the production of a transversely 
polarized p^ (including the interference from amplitudes 
with positive and negative helicity) by W^jyi while 
WijY results from the interference between longitudinal 
and transverse p^ polarizations. 

The terms W^^ can be expanded (sec Eqs. 4.10 and 
4.17 of Ref. [ll[) into trigonometric functions of the an- 
gles (t)s,<P, and (/?, where the coefficients are SDMEs (or 
combinations thereof) for W^^ , and n^^, and 



^ Because the target polarization is transverse to the incoming 
beam, there is a small longitudinal polarization with respect to 
the direction of the virtual photon. The effect of the latter and 
of a small longitudinal polarization of the beam be discussed 
later. 



for W^^ . Here the letters u,n, and s stand for unpolar- 
ized, normal, and sideways (with respect to the direction 
of the virtual photon and the electron scattering plane) 
target polarization, and the sub(super)scripts refer to the 
helicity of the virtual photon {p^ meson) in the helicity 
amplitudes that occur in the SDMEs. In the case of Wuu 
there are 15 independent terms in the expansion. There 
is a direct relation between these SDMEs and the ones in 
the Schilling- Wolf formalism Q . For Wut the expansion 
contains 30 independent terms. 

First the 15 'unpolarized' SDMEs of Wuu were deter- 
mined by fitting the angular distributions of the com- 
bined events for the two target polarization states. The 
fit was performed by maximum-likelihood estimation 
with a probability density function 

}u[c^,V.^) = ^u^ A{cj,,^,^)Wuu[c^,V.^). (7) 

where the function A represents the acceptance of the 
Hermes spectrometer. The factor Mu represents the 
normalization integral of the probability density func- 
tion, which was computed numerically using Monte Carlo 
events that are within the acceptance of the spectrome- 
ter. The non-exclusive background was included in the 
fit function using fixed effective values of the SDMEs for 
this background. The latter were obtained from a fit of 
the angular distribution of the Pythia6 Monte Carlo 
events for Ai? < 0.6 GeV. The results for the 15 unpo- 
larized SDMEs, which as mentioned arc for data taken in 
the years 2002 — 2005, arc fully consistent with those from 
the analysis of all data taken in the period 1996 — 2005 
using the Schilling- Wolf formalism [l^] . 

Then the 30 SDMEs Wut were determined, keeping 
the unpolarized SDMEs fixed to the values found in the 
fit of Wuu described above, using the probability density 
function 

/t(0,0s,</',i9) - AAj;i^(0,0s,<^,^) X 

{Wuu{^, V, ^) + PtWut (A 0s, ^)) ■ 

(8) 

A factor dip / d(j)s , which takes into account that the yields 
are evaluated differentially in the angle (ps, rather than 
in ■0, was left out, since it was very close to unity. As 
in the unpolarized case, the background was included in 
the fit. Since nothing is known about the asymmetry 
of the background, the 30 SDMEs for the background 
were taken to be zero, and the possible influence of this 
assumption was included in the systematic uncertainties. 

Besides the target polarization, various other sources 
of systematic uncertainties for the SDMEs and asymme- 
tries extracted were investigated and evaluated. In most 
cases the resulting systematic uncertainties were found 
to be negligible, i.e., very small compared to the statis- 
tical uncertainty. These include the effect of radiative 
corrections and the uncertainties resulting from the un- 
certainty in the unpolarized SDMEs and the background 
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FIG. 4: Values of SDMEs, or combinations thereof, for a transversely polarized proton target and an unpolarized beam. The 
SDMEs are sorted into three categories, which are separated from each other by the solid horizontal lines. From top to bottom: 
SDMEs containing s-channel helicity- conserving amplitudes, combinations containing at least one s-channel helicity-changing 
amplitude, and SDMEs containing two s-channel helicity-changing amplitudes. Within the second category the combinations 
are sorted into three groups associated with different virtual photon and polarizations. The inner error bars represent the 
statistical uncertainties. The full error bars represent the quadratic sum of the statistical and systematic uncertainties. In 
addition there is an overall scale uncertainty of 8.1% due to the uncertainty in the target polarization. 
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fraction. The uncertainty due to the angular dependence 
and asymmetry of the background was taken as the dif- 
ference between a fit with a background with no angular 
dependence and asymmetry, and one having the same 
angular dependence and asymmetry as the data. The 
resulting uncertainty was found to be negligible. 

The influence of the net beam polarization of approxi- 
mately 0.095 was estimated by including the SDMEs for 
Wlu and Wlt in the fit. Even if the latter had large 
uncertainties, the infiuence on the ones for Wut was 
negligible. The data presented in Fig. [5] are effectively 
integrated over all or two of the variables Q^, xb, and 
t' within the experimental acceptance. The effect of this 
kinematic averaging was estimated by comparing the re- 
sults of a Monte Carlo simulation that included a mod- 
elled dependence of the asymmetry on these variables 
with the model input values at the average kinematics. 
Also this effect was found to be negligible. 

In the extraction of the SDMEs the small longitudi- 
nal component of the target polarization with respect to 
the direction of the virtual photon (the average value 
of \Sl/Pt\ was 0.072) was neglected. This component 
introduces a term SlWjjl, which is described by 14 
SDMEs. As the value of Sl is small, these SDMEs can- 
not be determined from the present data. A system- 
atic uncertainty was estimated by using several sets of 
random values obeying the positivity bounds given in 
Ref. for these SDMEs, and evaluating the resultant 
changes. Changes of on average 55% of the statistical 
uncertainty were found, with a maximum of 76% for one 
SDME (Im(.s^^ + esQ (5")). This is the main source of 
systematic uncertainty. 

Lastly there are systematic uncertainties arising from 
misalignment of the detector, detector smearing effects, 
and bending of the beam and produced charged particles 
in the transverse holding field of the target magnet. The 
uncertainties due to all effects together were investigated 
with a Monte Carlo simulation of the possible influence 
of these effects. The resultant uncertainty was found to 
be negligible. 

The resulting SDMEs are shown in Fig. 21 Almost 
all of them are consistent with zero within 1.5(t, where 
tr represents the total uncertainty in the value of an 
SDME. Note that these include s-channel helicity con- 
serving SDMEs. Similar SDMEs in the unpolarized case 
were found to be non-zero and large (0.4 - 0.5). 

The SDMEs lm{s^+ - So + )' ImsZ|, and lmn^° de- 
viate more than 2.5a from zero. The former two involve 
the interference between natural (N) and unnatural (U) 
parity exchange amplitudes For instance, Ims"^ 

contains the product Nq^{U^^)* and ImsZ^ contains 
the product NZ^{U^^)* . The detailed analysis of unpo- 
larized data has shown that N^^ and NZ^ are dominant 
N amplitudes. The U amplitudes presumably are small, 
as they are suppressed at large Q^. However, is rel- 
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FIG. 5: The extracted amplitudes of the sm{(j) — (ps) com- 
ponent of AuT for longitudinally (top) and transversely po- 
larized (bottom) mesons. The inner error bars represent 
the statistical uncertainties. The full error bars represent 
the quadratic sum of the statistical and systematic uncer- 
tainties. In addition there is an overall scale uncertainty of 
8.1% from the uncertainty in the target polarization. 



atively large PJ,|2fl|. The SDME Imn^^ corresponds to 



a 7^ Pl transition, the SDMEs of which were found 
to be non-zero in the unpolarized case. The value of 
— 0.069 ± 0.022 measured for ImriQ ° is another indica- 
tion of violation of SCHC in the "fx ^ Pl transition. 

As mentioned, the sin((/) — (f>s) term in the transverse 
target-spin asymmetry for production of longitudinally 
polarized p° mesons is of special importance because of 
its sensitivity to the GPD E. The amplitude of this term 
is given in terms of SDMEs as [lll | 



,LL,E 



a{<p-<ps) 



Im(nOO 



- en 



, 



, 



, 
'0 



(9) 



The resultant values for all selected data and for bins in 
X, Q^, and t' are shown in Fig. [5] (top). They are all 
zero within the error bars. Because the SCHC violat- 
ing terms Im(n^^ ) and in Eq. [9] require a double 



helicity flip (see Ref. [Uj for details), they presumably 



LL,sin(</)— 05 ) 
UT 



can be neglected. Then the value of A 
—0.035 ± 0.103 ^ can be compared to the results of GPD 
calculations for the production of a longitudinally po- 



This is the value for 'all' data, which has average kinematics 
< q2 >= 1.95 GeV^, < Xb >= 0.08, and < -t >'= 0.13 GeV^. 
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which is 



larizcd by a longitudinal photon A^^t^, '^^^ 
given by lTa{n^^)/u^^ . 

The sin(0 — (ps) amplitude for the production of trans- 
versely polarized p° mesons is given by 



A 



TT,sin{<p-<j>s) 
UT 



Im(?' 
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'0 , 



1 - (u 





++ 



■ eu 



(10) 



, 



The values for this asymmetry are shown in Fig. [5] (bot- 
tom). Also these are zero within the error bars. 

A few groups have performed GPD-based calcula- 
tions of the transverse target asymmetry for exclusive 
production. In Refs. [1 [2l| the quark GPD S"? is 



taking J'^ = 0. 



P 

parametrized in terms of the value of 
Ref. 2l| includes the contribution of gluons. The calcu- 
lated values oi Af^rl^t,~:t'} 

0, 



^UT fi PL range 0.15 to 0.00 for 

J" = 0.0 to 0.4. In Refs. [H, HI the GPDs are modeled 



using data for nucleon form factors, sum rules and pos- 
itivity constraints. The results of both calculations are 
similar. Values of J" and J'^ of approximately 0.22 and 
0.00, respectively, are found, and the calculated values 
of the asymmetry are very small (—0.03 to 0.02), which 
is consistent with the present data. It must be realized 
that the results depend on the modelling of the relevant 
GPDs of both quarks and gluons, and that the kinematic 
conditions of the calculations are in several cases outside 
the kinematic range of the present data. 

In summary, the transverse target single-spin asymme- 
try was measured for exclusive electroproduction on 
a transversely polarized hydrogen target. Spin density 
matrix elements were determined by using the angular 
distributions of the produced p° mesons and their de- 
cay into two pions. Almost all of the SDMEs describing 
transverse target polarization were found to be consis- 
tent with zero. A notable exception is an SDME that 
corresponds to the production of a longitudinally polar- 
ized p° by a transverse virtual photon. The fact that it 
is non-zero indicates a small violation of s-channel helic- 
ity conservation in the case of a transversely polarized 
target. The amplitude of the sin((/) — (j)s) component of 
the asymmetry for the production of longitudinally polar- 
ized p° mesons was found to be small (—0.035 ± 0.103). 
Neglecting double helicity changing SDMEs, this com- 
ponent can be identified with the leading-twist term of 
the asymmetry. Calculations based on generalized par- 
ton distributions predict small values, consistent with the 
measured value. 
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